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We discuss the status of different determinations of as, motivating a precise and reliable computa- 
tion from lattice QCD. In order to suppress perturbative errors, the non-perturbative computation 
has to reach high energy scales /i. Such results already exist in the SF-scheme for A^f = 0,2 [1,2] 
and Nf = 3 [3]. We recently added the running with four massiess flavours in a range of a from 
about 0.07 to 0.3 . It is based on our recent determination of the Sheikholeslami Wohlert coeffi- 
cient in the four-flavour theory. 
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1. Introduction 

The strong coupling, Us, represents a fundamen- 
tal parameter of the strong interactions. Its scale de- 
pendence in a suitable renormalization scheme teaches 
us about the connection between the strongly non- 
perturbative and the dominantly perturbative regions 
of QCD. Its value at high energy is needed for phe- 
nomenology, for example for the prediction of Higgs 
production cross sections for the LHC. The uncertainty 
on tts is not negligible in this context [6] . 

For this conference we had compiled a plot of pre- 
cise determinations of a = ttg with the renormaliza- 
tion scale set to the mass of the Z-boson and in the 
MS scheme. It is shown in Fig. 1. We do not claim 
completeness but rather want to illustrate that the many 
determinations do not agree well within the estimated 
uncertainties. 

Indeed, the spread of results in Fig. 1 is not that 
surprising, since a precision determination of as{Mz) 
is difficult and so far compromises on various sources of eiTors had to be made despite an ever 
increasing sophistication in the analysis. Sources of errors are: 

• Low energy: most determinations are not done from a process with an energy scale of order 
Mz but at significantly lower energies and are then evolved perturbatively to jj. = Mz- A 
prominent example is the determination from T-decays, labeled R^;. 

• Sophisticated analysis involving simultaneous fits to OCs(At) and non-perturbative parameters 
of QCD, such as structure functions [6] and parameters of SCET [7]. 

• Global fits to many processes (EW fits) [8], which of course means that the correctness and 
theoretical mastering of the standard model enters in detail. 

• The use of bare, unrenormalized, perturbation theory of Wilson loops at the cutoff scale in 
the lattice determination [4] shown in the figure. We discuss this further below. 

• None of these applies to the extraction from the hadronic cross section on the Z-peak, Rz, 
but this is experimentally more difficult, resulting in a larger quoted error. 

The present world average by S. Bethke is dominated by the lattice determination [4]. We find 
this worrying due to several reasons. First there is the use of bare perturbation theory. There the 
problem is that successive terms do not show a convincing 'convergence' pattern. Either there are 
large coefficients of higher order terms (expansion in the original bare coupling) or the (modified) 
coupling itself "lives" at a quite low energy scale where it is lai^ge ("tadpole improved" coupling, 
coupling in "potential scheme"). In either case errors due to left out remainder terms are difficult 
to quantify. The continuum limit is in this way only reached at an asymptotic rate proportional to 
I/I lna| and can in practise not be taken. Furthermore, rooted staggered quarks with their doubtful 
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Figure 1: Some recent precise numbers 
for as(Mz) in the MS scheme. From top 
to bottom: lattice QCD by HPQCD [4], 
hadronic T decays (7?t) [5], deep inelas- 
tic scattering [6], thrust [7], global elec- 
troweak fit [8], hadronic decays of the Z 
(Rz) [9, 10] and the world average of [1 1]. 
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theoretical status are employed. In spite of this, a phenomenological observation adds support to 
the analysis of [4] : there is agreement of results from several different variations of the theme. 

The above mentioned difficulties have been a motivation for the ALPHA collaboration to work 
on a program which essentially is based only on the assumptions that the continuum limit of the 
lattice theory exists and asymptotic freedom is present non-perturbatively. The resulting errors in 
tts are dominantly statistical and also the systematic component can be reduced further when the 
overall precision is improved in the future. In particular, as has been explained many times [12-14], 

• a is defined non-perturbatively in a physical (regularization independent) scheme, 

• this running can be efficiently computed numerically by a finite size strategy, 

• the continuum limit can be taken in individual steps free of multi-scale problems, 

• renormalized perturbation theory is used only at large scales, where its precision is further- 
more validated non-perturbatively. 

These properties come with a price. Since a specific scheme had to be devised in order to make 
precise non-perturbative computations possible, the j3 -function at 3-loop order needed to be cal- 
culated [15-22]. The same perturbative calculations are also important to remove dominant parts 
of the discretization errors (see section II.2 of [14]). Furthermore a controlled non-perturbative 
running clearly has to include the strange and charm quarks; this step has only been carried out 
recently and we report on it here. Still, as will be discussed in the conclusions more work is needed 
before a numerical value for as{Mz) can be given for phenomenology. It is also for these reasons 
that the ALPHA collaboration has so far not been able to publish a value of a in the physical theory 
with all quarks which we would like to be taken into account in the world average. 

2. Computation of the step scaling function 

In our finite size strategy, the coupling 
is defined in a Euclidean space-time of 
size with Schrodinger functional bound- 
ary conditions [15,16] and a renormalization 
scale jx = l/L. The discrete scale evolution 
defines the step scaling function o via 

g'{sL) = a{s,f{L))= lim Z{s,g'{L),a/L) 

a/L-^Q 

which can be computed non-perturbatively 
as the continuum limit of the lattice approxi- 
mant Z as indicated. We chose a scale factor 
5 = 2 for practical reasons and omit this ar- 
gument from now on. In the Schrodinger functional, Dirichlet boundary conditions are imposed 
in Euclidean time and therefore the 0(a) Symanzik improvement of the theory requires boundary 
terms. Their coefficients are taken from perturbation theory, exactly as in [23], where we deter- 
mined the crucial bulk 0(a) improvement coefficient Csw non-perturbatively. 
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Figure 2: Interpolation eq. (2.1). 
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The computation of the lattice step scaUng function, l.{u,a/L) requires to tune the bare mass 
and the bare couphng of the theory such that the PCAC mass (defined exactly as in [2,24]) van- 
ishes and g^{L) = u. At the same bare parameters one then computes 'L{u,a/L) = g^{2L). An 
explicit 2-dimensionaI tuning is rather cumbersome. We therefore followed [25], picked a se- 
ries of bare couplings = 6/(5 and tuned the PCAC mass to zero. For those bare parame- 
ters we then compute g^{p,L/a) and g^{p,2L/a) and interpolated to the desired values of u via 



gHP,L/a) 



M 
m=0 



-1 



3.5 



(2.1) 

motivated by perturbation theory. We do not 
fix the known perturbative expansion coeffi- 
cients, not even Co^/a- The interpolations for 
L/a = 8 are illustrated in Fig. 2. £(«, 1/8) 
is given by the value of the upper curve at 
the j3 where the lower one passes through 
g^(j8,8) = u. Stability of the interpolations 
with respect to M was checked. 

Using £ from the interpolation, we form 
the 2-loop improved lattice step scaling func- 
tion [26] 
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with 81,82 known from [19-22,27]. We ex- 
pect £(^) to have smaller overall cutoff ef- 
fects. Asymptotically, they still start at or- 
der a X but terms of order a'" x u" are 
removed for all m and for n < 3 (in fact 
non-perturbatively in a). As mentioned pre- 
viously, the order a x terms are due to 
the only perturbatively known boundary im- 
provement terms. Their influence was ex- 
plicitly checked for Nf = 2 and found to be minor [2, 14] for our action, at least when Ct is known 
to 2-loop order. We therefore assume that the step scaling function converges effectively at a rate 



Figure 3: Constant fit continuum extrapolation. 



(2.3) 



To study the continuum limit and its uncertainty we caiTied out three different analysis. 

• Constant jit: K^tofTP'\u,a/L) for L/a = 6, 8 to a constant, for each m. 

• Global fit: A fit L^^\u,a/L) = a{u) + pu'^ i^/L')^ : with a separate, independent parameter 
a{u) for each value u but a common parameter p modelling the cutoff-effects. 

• L/a = 8 data: Using directly o{u) = L^^\u, 1/8). 
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The three different ansatze yield results which are in complete agreement with each other. The 
global fit returns p = 0.007(85) which is a good indication that cutoff effects are negligible in the 
data for L/c? = 6, 8. In order to have a safe error estimate on the continuum hmit we chose just the 
L/a = 8 data as our present result. 

3. The running of the coupUng 

A polynomial interpolation a{u) = u + squ^ + siu^ +0.0036m^ — 0.0005 m^, with the coeffi- 
cients up to fixed by perturbation theory represents a{u) in the range 0.9 <u< 2.7 with negli- 
gible interpolation errors. The running of the coupling is then obtained by solving the recurrence 

Ui = a{ui+i), i = 0,...,n, Uo = Umax = f (Lmax) (3.1) 

for Ui = g^{Li), Li = 2~'Liiiax- As Fig. 4 shows, agreement with perturbative running at the 3-loop 
level is found at the highest scales in agreement with standard estimates of remainder terms. This 
allows us to relate to L9 x A by using the 3-loop /3 -function [for couplings up to M9 only]. Then 
with the non-perturbative g{u) used in eq. (3.1) we can connect to larger values of u, for example 

ln(AL„,ax) = -2.294(83) at f{Lm,x) = "max = 3.45 . (3.2) 

While a precise MeV value still has to be determined we clearly expect L^ax to lie in the range 
of hadronic scales. The (uncoiTclated) errors of our primary MC data for g^(j3,L/a) are propa- 
gated through all steps of the analysis as described in [24]. Note that some steps, in paiticulai^ 
eq. (3.1), introduce correlations into the final results. The eiTors of the points Fig. 4 are therefore 
not independent. 

4. Conclusions 

At the present level of sta- 
tistical errors, systematic errors 
from the continuum extrapola- 
tion are almost certainly negli- 
gible. In the future, in partic- 
ular with further reduced sta- 
tistical errors, data with larger 
L/a clearly have to be added. 
It will also be interesting to 
compare the efficiency of com- 
putations with different regu- 
larizations of the Schrodinger 
functional, such as chirally ro- 
tated boundary conditions [28, 
29] and staggered quarks [30, 
31]. 
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Figure 4: Non-perturbative running coupling compared to pertur- 
bation theory. 



5 



Running of the SF-coupling with four massless flavours 



Rainer Sommer 



We do observe a small but significant deviation from 3-loop perturbation theory at the largest 
coupling reached in Fig. 4. It is about 10% (three standard deviations) and the Schrodinger func- 
tional coupling has a value of ttsF ~ 0.28. For Nf = 2 a. similar effect was visible only for larger 
coupling [2]. These findings underline the necessity of going to weak coupling before applying 
perturbation theory. 

We are now a good step closer to the computation of the A-parameter in 4-flavor QCD, which 
may then be perturbatively connected to e.g. the 5-flavor aj^(Mz). However the low energy 
scale Lniax that was introduced for technical convenience, remains to be expressed in physical units 
through large volume 4-flavor simulations and we may want to improve the precision in Fig. 4. 
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